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Abstract 

With the current burst of network theory research (especially in connection 
with social and biological networks) there is a renewed interest on realiza- 
tions of given degree sequences and uniform sampling of those realizations. 
In this paper we propose a new degree sequence problem: we want to find 
graphical realizations of a given degree sequence on labeled vertices, where 
certain would-be edges are forbidden. Then we want to sample uniformly 
all possible realizations. 

In this paper, as a first step, we solve this restricted degree sequence 
(or RDS for short) problem if the forbidden edges form a bipartite graph, 
which consist of the union of a (not necessarily maximal) 1-factor and a 
(possible empty) star. Then we show how one can sample the space of 
all realizations of those RDSs uniformly if the degree sequence describes 
half-regular bipartite graphs. Our result contains, as special cases, the well- 
known result of Kannan, Tetali and Vempala on sampling regular bipartite 
graphs and a recent result of Greenhill on sampling regular directed graphs 
(so it also provides new proofs of them). 

The RDS problem with one forbidden 1-factor and with one star is 
self-reducible, therefore our fully polynomial almost uniform sampler (a.k.a. 
FPAUS) on the space of all realizations also provides a fully polynomial ran- 
domized approximation scheme (a.k.a. FPRAS) for approximate counting 
of all realizations. 

Keywords: restricted degree sequences; rapidly mixing MCMC; Sinclair's 



*We acknowledge financial support from grant #FA9550-12- 1-0405 from the U.S. Air 
Force Office of Scientific Research (AFOSR) and the Defense Advanced Research Projects 
Agency (DARPA). 

^Partly supported by Hungarian NSF, under contract NK 78439 

^Partly supported by Hungarian NSF, under contract K77476 and NK105645 

^Partly supported by Hungarian NSF, under contract PD84297 

^Partly supported by Hungarian NSF, under contract NK 83726 



1 



multicommodity flow method; FPAUS; self-reducible problem; FPRAS 



1. Introduction 

In the last fifteen years, network theory has been undergoing an exponential 
grow. One of its more important problems is to algorithmically construct 
networks with predefined parameters or uniformly sampling all networks 
with these given parameters. For general background, the interested reader 
can turn to the now-classic book of Newman, Barabasi and Watts ( [14] ) or 
to the more recent book of Newman ([15]). 

In this paper we study networks (or graphs as we like to refer them) with 
given degree sequences. We propose a new degree sequence problem class 
called restricted degree sequence problem and solve its first instance: 
we build procedure to decide quickly whether there exists a graph G with 
a given degree sequence, where G completely avoids a predefined set of for- 
bidden edges, then develop a fast mixing Markov Chain approach to sample 
almost uniformly all different realizations of the sequence. 

We will focus on the technical details of this problem; we do not intend 
to survey the history of the degree sequence problems in general or their 
connections to other developments in network theory. One can find detailed 
information about this specific background in the recent paper of Greenhill 
([5]), or in our previous paper [13]. Therefore we will touch only the directly 
affected definitions and earlier results. 

In the next section we discuss some of the known degree sequence prob- 
lems and introduce our proposed new problem class. In Section 3 we will 
study in details one instance of this new problem class, describing and prov- 
ing a fast Havel-Hakimi type greedy algorithm to construct realizations. 
Then in Section 4 we discuss some known MCMC approaches to sample 
different kind of degree sequence realizations and state our result about the 
rapid uniform sampling of this instance of our proposed model. In Sections 
5 and 6 we will describe our approach in details. In Section 7 we show that 
the studied problem is a self-reducible one, therefore our almost uniform 
sampling method provides good approximate counting of all realizations. 
Finally in the Appendix we will discuss a necessary technical detail of the 
required, but very slight generalization of the "Simplified Sinclair's method" 
introduced in [13]. 

2. Degree sequences and restricted degree sequences 

Let's fix a labeled underlying vertex set F of n elements. The degree sequence 
d(G) of a simple graph G = iV^E) is the sequence of its vertex degrees: 
d(G)i = d{vi). (Here "simple" means that there are no loops or multiple 

edges. In general, in cases where multiple edges and/or loops are not ex- 
cluded the corresponding results are easier. See for example Ryser, [17].) A 
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non-negative integer sequence d = {di, . . . , (i„) is graphical iff d(G) = d for 
some simple graph G, and then G is a graphical realization of d. The simplest 
method (a straightforward greedy algorithm) to decide the graphicality of 
an integer sequence was discovered by Havel ([7]), and was rediscovered in- 
dependently later by Hakimi ([6]). The validity of their algorithm is proved 
by the means of the swap operation. It is defined as follows: 

Let G be a simple graph and assume that a,b,c and d are different 
vertices. Furthermore, assume that (a, c),(6, d) G E{G) while {b,c),{a,d) 
E{G). Then 

EiG') = E{G) \ {(a, c), (6, d)} U {(6, c), (a, d)} (2.1) 

is another realization of the same degree sequence. We call this operation 
as swap and will denote by ac, hd =^ be, ad. 

The analogous notions for bipartite graphs are the following: if i? is a 
simple bipartite graph then its vertex classes will be denoted by U{B) = 
{ui, . . . , tifc} and W{B) = {wi, . . . , w^}, and we keep the notation V{B) = 
U{B) U W{B). The bipartite degree sequence of B, hd{B) is defined as 
follows: 

hd{B) = ((d(txi), . . . , d{uk)) , {d{wi), d{we))y 

We can define the swap operation for bipartite realizations similarly to (2.1) 
but we must take some care: it is not enough to assume that (6, c), (a, d) ^ 
E{G) but we have to know that a and b are in the same vertex class, and c 
and d are in the other. 

To make clear whether a vertex pair can form an edge in a realization 
or not we will call a vertex pair a chord if it can hold an actual edge 
in a realization. Those pairs which cannot accommodate an edge are the 
non-chords. (Pairs from the same vertex class of a bipartite graph are 
non-chords.) 

For directed graphs we consider the following definitions: Let G denote a 
directed graph (no parallel edges, no loops, but oppositely directed edges 
between two vertices are allowed) with vertex set X{G) = {a;i,X2, • • • ,iCn} 
and edge set E{G). We use the bi-sequence 

dd(G) = ( 4, . . . , 4) , (dr, a!2-, . . . , ) ) 

to denote the degree sequence, where df denotes the out-degree of vertex 
Xi while d^ denotes its in-degree. A bi-sequence of non-negative integers is 
called a directed degree sequence if there exists a directed graph G such that 
(d+,d~) = dd(G). In this case we say that G realizes our directed degree 
sequence. 

We will apply the following representation of the directed graph G : let 
B{G) = (U,W;E) be a bipartite graph where each class consists of one 
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copy of every vertex. The edges adjacent to a vertex Ux in class U represent 
the out-cdgcs from x, while the edges adjacent to a vertex iv^ in class W 
represent the in-edges to x (so a directed edge xy is identified with the edge 
UxWy). If a vertex has no in- (out-) degree in the directed version, then we 
delete the corresponding vertex from B{G). (This representation is an old 
trick, applied already by Gale [4].) There is no loop in om directed graph, 
therefore there is no {ux,Vx) type edge in its bipartite realization - these 
vertex pairs are non-chords. 

In this paper we will study the following common generalization of all the 
previously mentioned degree sequence problems: 

The restricted degree sequence problem d"^ consists of a degree sequence 
d and a set C (^) of forbidden edges. The problem, as it is in the origi- 
nal problem, is to decide wether there is a simple graph G on V, completely 
avoiding the elements of J^, which provides the given degree sequence, fur- 
thermore to design a uniform sampler of all realizations. 

The bipartite restricted degree sequence problem bd^^ consists of a 
bipartite degree sequence bd on (^7, W), and a set J" C [U, W] of forbidden 
edges. The problem, as it is in the original problem, is to decide wether 
there is a bipartite graph G on {U, W) completely avoiding the elements of 
T, which provides the given degree sequence. 

Clearly a bipartite restricted degree sequence problem bd"^ on (U, W) is 
a restricted degree sequence problem d-^ on J7 U W, where J-"' = F U [J7]^ U 

[w?. 

It is easy to see that the restricted degree sequence problem is very 
closely related to Tutte's /-factor theorem ([19]). However, while Tutte's 
approach provides a polynomial time algorithm to decide wether a given 
degree sequence satisfies a d-^ problem, it does not provide an easy greedy 
algorithm to do so, furthermore it does not help in the sampling problem. 

As we already mentioned the RDS notion is a common generalization 
of several "classical" degree sequence problems. For instance, the case of 
bipartite degree sequences is such an example. Another well known example 
is the case of directed degree sequences: in its bipartite representation, a 
forbidden 1-factor excludes the directed loops in the original directed graph. 

Definition 2.1. Let d"^ be a restricted degree sequence problem and let 

G be a realization of it. The sequence of vertices C = {xi,X2, ■ ■ ■ ,X2i) is 
a chord-circuit if: (Dl) all pair xiX2,X2X3, . . . X2i-iX2i,X2iXi are chords; 
and (D2) each chord occurs only once. 

A chord-circuit is elementary if (D3) no vertex occurs more than twice; 
furthermore (D4) when two copies of the same vertex exist, then their dis- 
tance along the circuit is odd. 

The chord-circuit C alternates if the chords along C are edges and non- 
edges in turn (for example X2j-iX2j are edges for 1 < j < i, while the other 
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chords are not edges in G). 

Deleting the actual edges along C from G and adding the other chords as 
edges constructs a new graph G' which is again a realization of d-^. This 
operation is known as a circular C2i-swap and denoted by Sc- 
Finally, two different vertices of the alternating chord-circuit C form a po- 
tential vertex pair (or PV-pair for short) if (i) this chord docs not belong 
to the circuit; and (ii) the distance of the vertices distance along the sequence 
(which is the number of chords between them) is odd and it is not 1. If each 
PV-pair is a non-chord (that is G T), then this circular swap is called a 
J^-compatible swap or T-swap for short. 

The T-swap is one of the central notions of this paper. When i = 2 then 
the circular C4-swap coincides with the classical Havel-Hakimi swap. When 
i = 3 then we get back the notion of the triangular Cg-swap, which was 
introduced in paper [2] in connection with directed degree sequences. 

We define the weight of the J^-compatible circular C2i-swap as w{C2i) = 
i — 1. This definition is in accordance with the definitions of the weight of 
the classical HH-swaps, and the weight of a triangular Ce-swap in paper 
[2]. Furthermore it is well known (see for example again [2]) that in case of 
simple graphs (i — 1) Havel-Hakimi swaps are needed to alternate the edges 
along C2i ■ As we will see next the same applies for any (elementary) circular 
C2i-swap: 

Lemma 2.2. Let G be a realization of d-^ and let the elementary chord- 
circuit C of length 2i alternate. Then the circular C-swap operation can be 
carried out by a sequence of -swaps of total weight i — 1. 

More precisely there exists a sequence G = Go,Gi, . . . ,G£ of realizations 
such that for each j = 0, . . . ,£ — 1 there exists an J^-swap operation from Gj 
to Gj+i, the difference between G and Ge is exactly the alternating circuit 
C, finally the total weights of those J^-swap operations is i — 1. 

Proof. We apply mathematical induction for the length of the chord-circuit: 

assume this is true for all circuits of length at most 2i — 2. (For i = 3 we 
already saw this.) Then take an alternating elementary chord-circuit C of 
length 2i in a realization of d-^. 

If each PV-pair in C is a non-chord, then the circular C2i-swap itself is 
a J^-swap of weight i — 1. So we may assume that there is a PV-pair uv in 
C which is a chord. This chord together with the two "half-circuits" of C 
form chord-circuits Ci and C2 using the chords of the original circuit C and 
twice the chord uv. One of them, say Ci, is alternating. The length of Ci is 
2j < 2i therefore there exists a J^-compatible swap sequence of total weight 
J — 1 to process it. After the procedure the status of uv will alter into the 
other status. With this new status of the chord the circuit C2 becomes an 
alternating one, so it can be processed with _ 1 total weight - and 
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after this procedure the chord uv is switched back to its original status. So 
we found a swap sequence of total weight i — 1 which finishes the proof. □ 

The space of all realizations: Consider now the set of all possible real- 
izations of a restricted graphical degree sequence d-^. Let G and H be two 
different realizations. The natural question is whether G can be transformed 
into H using J^-swaps? 

For classical degree sequences this problem was solved affirmatively al- 
ready in 1891 by Petersen in the (by now almost completely forgotten) 
paper [16]. Havel's paper [7] also provided (an implicit) solution. For bi- 
partite graphs (with possible multiple edges but no loops) this was done by 
Ryscr ([17]). For simple bipartite graphs it was folklore. Finally for directed 
graphs it was done in [11] (and later rediscovered in [3]). 

Theorem 2.3. The space of all realizations G = (V, E) of the restricted 
degree sequences problem d-^ is connected. Therefore the usual Markov chain 
defined on G is irreducible. 

Proof. What we have to prove is the following: let G and H be two realiza- 
tions of d"^. Then we have to find a series of realizations G = Gq, . . . , Gi-i, 
Gi = H, such that for each j = 0, . . . , i — 1 there exists an J-'-swap from Gj 
to Gj+i. (Such realization pairs form the edge set of the Markov chain G.) 

Consider the symmetric difference of the edge sets of the two realizations: 
A = E{G)AE{H). This set is two-colored by the original hosts of the edges: 
there are G-edges and if-cdges. It is clear that in the graph Q = (V, A) 
for each vertex v the numbers of G-edges and i7-edges incident to v are 
the same: dciv) = duiy). By Euler's method it can be decomposed into 
alternating circuits Ci, . . . , Q. 

Let's recall that a circuit in a simple graph G is a sequence of vertices 
^0; • • • ! V2t, where ?7o = ^2* s.t. the consecutive vertices are adjacent and each 
edge can be used at most once. Note that there can also be other indices 
i < j such that vi = vj. A circuit is called a cycle, if it is simple, i.e., for 
any i < j, Vi = Vj only if i = and j = 2t. A circuit is alternating for G 
and H if the edges come in turns from E{G) and E{H). When this is the 
case then the corresponding chord-circuit in realization G (as well in H) is 
alternating. 

We can find a decomposition that no circuit contains a vertex v twice 
s.t. their distance 6 (the number of edges between them) is even. Indeed, 
since Q is simple, therefore 6 is at least four, consequently the vertex v can 
split the original circuit into two smaller, but still alternating circuits. So 
for example an alternating circuit decomposition with maximal number of 
circuits has this property. It also implies that no circuit may contain a 
vertex three times otherwise at least two copies of the vertex would be of 
even distance from each other. 



6 



The application of Lemma 2.2 proves that each circuit C can be processed 
with |C|/2 — 1 total weight. This finishes the proof. □ 

In the paper [2] the following formula was developed for the required 
minimum weight of transforming one realization in to an other one of the 
same (unrestricted) degree sequence: Consider the realizations G and H and 
the symmetric difference EAF of their edge sets. Denote maxC(G, H) the 
maximum possible number of the alternating circuits in such a decomposi- 
tion. 

Theorem 2.4 (Erdos - Kiraly - Miklos [2]). For any pair G, H of realizations 
of a given degree sequence the weight of the shortest swap sequence between 
the two realizations is 

dist(G, H) = \EiG)A^m _ ^^^^^^ (2.2) 

Now we show that the analogous result holds for any restricted degree se- 
quence. 

Theorem 2.5. Let G, H be two realizations of the same restricted degree 
sequence d^. Then 

dist^(G, H) = \E(G)AE{H)\ _ ^axC(G, H). (2.3) 
Proof. We show that 

RHS =^ dist{G,H) <^ dist^(G,iJ) <^ RHS. (2.4) 

The first equality is just Theorem 2.4. To check the second inequality assume 
that the J^-swap sequence 

<5C:,...,5c„ (2.5) 

transfers G to H. If the length of Cj is 2k j, then, applying Theorem 2.4 
again, we have that Sc^ can be obtained as a composition of kj — 1 many 
"standard" swaps Sc^ i , • ■ ■ , Se^ _j . So the "standard" swap sequence 

Sci,i, ■ ■ ■ ,Scij^_j^,Sc2,i, ■ ■ ■ ,Sc2,j2-i ■ • •'5'c„,i, ■ ■ ■,Sc„j„-i (2.6) 

transfers G to H, and the weight of this sequence is Y^^=ii^j ~ !)> which 
is exactly the weight of the J"-swap sequence Sci, ■ ■ ■ ,Scn from (2.5). Thus 
the second inequality holds. 

Finally the third inequality holds by Lemma 2.2. □ 
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3. The star+factor restricted degree sequences 

We turn our interest now for the following specialized restricted degree se- 
quence problem: d-^ is called a Star — 1-Factor Restricted Degree 
Sequence problem (or star+factor problem for short), if 

{^) the set T of the forbidden edges is a bipartite graph where the edges 
are the union of an 1-factor and a star with center s. 

Similarly, if bd is a bipartite degree sequence, and (^') holds for T, then 
bd"^ is called a Bipartite Star — 1-Factor Restricted Degree Sequence 
problem (or bipartite star+factor problem for short). 

As we already mentioned, Tutte's /-factor theorem can always be utilized 
to find actual graphical realizations of our star+factor restricted degree se- 
quence. However in this special case we can apply a Havel-Hakimi type 
greedy algorithm to construct such realizations. 

Consider our star+factor degree sequence problem d-^. For a given ver- 
tex X & V denote C{x) the set of those vertices in V which form chords 
together with x. If 

(i) for each y ^ C{x) has at most one vertex, denoted by , such that 
chord yy-^ belongs to furthermore if y^ = then y = z, 

then we say that C{x) is normal, and we fix an order -<x on C{x) such that 

(ii) if d{y) > d{z), or d{y) = d{z) and d{y-^) > d{z-^), then y -<x z. 

Lemma 3.1. Let G be a graphical realization of our star+factor RDS , 
let X ^ V and assume that C{x) is normal with the order -<x- Assume 

furthermore that y z and xz £ E while xy ^ E. Then there exists an 
alternating chord-circuit C = (y, x, z,vi, . . . ,Vi) with i = 1 or 3. In the 
acquired new realization we have Tg'{x) = Tg{x) \ {z} U {y}. 

This statement is actually almost the same as Lemma 4 in [3] and its proof 
could be recalled. However here we give a complete proof. On one hand this 
keeps this paper self-contained, on the other hand paper [3] uses a different 
language. 

Proof. We have xz e E but xy ^ E. At first assume that there exists a 

vertex u ^ x,y, z, such that uy G E, and uz ^ E but u ^ z-^ . When such 
vertex exists then C = (x, z, ti, y) is a suitable alternating chord-circuit as 
the xz, uy =^ xy, uz swap shows. 

Definition 3.2. From now on the notation xz, uy =^ xj/, uz always means 
that xz, uy are edges, xy, uz are non-edge chords, and we create a new 
realization with the indicated swap. 
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We continue the proof: when d{y) > d{z) then this happens automati- 
cahy since y belongs to at most one forbidden pair. However, if d{y) = d{z) 
then it can happen that z-^y € E and 

\/u ^ x,y, z,y^ , z''^ we have uy^E^^uz^E. (3-1) 

It is important to observe that in this case y^ z E, otherwise some u would 
not satisfy (3.1) (in order to keep d{y) = d{z)). 

So the only case when we do not find automatically an appropriate swap 
with X, y and z is when d{y) = d{z), yz-^ is an edge and zy-^ is a chord but 
not an edge. In this case, we can find a u ^ y, z such that y-^u € E but 
z-^u E since d{y^) > d{z^). 

Now C = {y,x, Zjy-^ ,u, z"^ ,y) is the required alternating chord circle. 
See the figure below. The three line types denotes the edges, the chords 
which are non-edges, finally the forbidden non-chords. □ 




Lemma 3.1 provides the following easy Havel-Hakimi type greedy algorithm 
to construct at least one graphical realization of our restricted degree se- 
quence problem. 

HH-algorithm for a star-l-factor degree sequence problem. (Recall that the 
center of the forbidden star is denoted with s. Since a star can be empty, 
we can assume that s is always defined): 

(Hi) take an ordering ~<s on C{s) (which is normal) and connect s to the 
first d{s) vertices (with respect to -<s) of C{s). Delete s and update 
the degrees of the used vertices accordingly. 

(-^2) take the remaining vertices one by one and repeat the process. 

Theorem 3.3 (Generalized HH-theorem for the star-l-factor RDS problem). 

The d"^ star+factor restricted degree sequence is graphical if and only if the 
previous greedy HH-algorithm provides a realization. 

Proof. As usual, it is enough recursively apply Lemma 3.1. Since we start 
our algorithm with the vertex s therefore C{s) satisfies condition (i) therefore 
it is normal. When s is deleted, then the remaining forbidden edges are from 
the original 1-factor, so the lemma applies automatically at each subsequent 
step. □ 



9 



In the case of a bipartite degree sequence d the situation is very similar: 
we define the normahty of any C{x) formally the same way as before. The 
definition of the order -<x is also the same as before. 

Here it is interesting to observe, that if x is in class U than C{x) is 
subset of class W, the vertices y-^ and z-^ belong again to class U, so u € W 
(and, finally, the forbidden edges define those vertices are elements of J-"). 
Furthermore, as in any bipartite degree sequence problem, if we determine 
all edges adjacent to the vertices in U then we automatically placed all edges 
adjacent to the vertices in W as well. 

Theorem 3.4. Lemma 3.1 and Theorem 3.3 apply for the vertices of class 
U in case of bipartite star+factor restricted degree sequences without any 
changes. 

Proof. Indeed, by definition, the center s of the forbidden star belongs to U . 
Also by definition, in the proof of the Lemma the vertices we are considering 
are from the vertex class U only: x € U. Therefore C{x) is a subset of W. 
Consequently the vertices y-^ and z-^ belong to vertex class U therefore the 
vertex u must belong to vertex class W again. Therefore all forbidden edges 
considered in the proofs belong to J^. So the proofs apply without changes 
for the bipartite case as well. □ 

4. Sampling uniformly half-regular, bipartite star+factor restricted 
degree sequences 

In this section we describe the main result of our paper which is an MCMC 
algorithm for (almost) uniform sampling of the space of all realizations of 
the 

Half-regular, Bipartite Star — 1-Factor Restricted Degree Sequence 

problem: Let d be a bipartite degree sequence with a star+ 1-factor type 
forbidden edge set T, where the center of the star is denoted by s (and 
belong to U). Furthermore let 

($) the degree sequence d is a half-regular bipartite one: it requires a 
B = (U, W; E) bipartite graph where each vertex u U - with the 
one possible exception s - has the same degree. We will write V for 
ULiW. 

We will apply a Markov Chain Monte Carlo method for almost uniform sam- 
pling of all possible realizations of our d-^. Originally the MCMC method for 
realizations' sampling was proposed by Kannan, Tetali and Vempala (1999, 
[9]). They conjectured that their process is rapidly mixing on the realiza- 
tions of any (unrestricted) degree sequences, i.e., starting from an arbitrary 
realization of the degree sequence, the process reaches a completely random 
realization in reasonable (i.e., polynomial) time. They managed to prove 
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the result for bipartite regular graphs. Their conjecture was proved for ar- 
bitrary regular graphs by Cooper, Dyer and Grccnhill (2007, [1]). (Their 
result does not automatically generalize the previous result, since their ver- 
sion does not allow forbidden edges.) An analogous theorem was proved 
by Grccnhill on regular directed graphs ([5]). Miklos, Erdos and Soukup 
proved in [13] that this Markov process is also rapidly mixing on each bi- 
partite half-regular degree sequence (here there is no exceptional vertex 
s). In this paper we will prove that the analogous Markov process is rapidly 
mixing on the half-regular star-l-factor bipartite restricted degree sequence 
problem. 

The state space of our Markov chain is the graph G = {V{G), E{G)) 
where V{G) consists of all possible realizations of our problem, while the 
edges represent the possible swap operations: two realizations (which will 
be indicated by upper case letters like X or Y) are connected if there is a 
valid J^-swap operation which transforms one realization into the other one 
(and the inverse swap transforms the second one into the first one as well). 

The transition (probability) matrix P is defined as follows: let the current 
realization be G. Then 

(a) with probability 1/2 we stay in the current state (that is, our Markov 
chain is lazy); 

(b) with probability 1/4 we choose uniformly two- two vertices ui,U2;vi,V2 
from classes U and W respectively and perform the swap if it is possible; 

(c) finally with probability 1/4 choose three - three vertices from U and W 

and check wether they form three pairs of forbidden chords. If this is 
the case then we perform a circular Ce-swap if it is possible. 

Here cases (b) and (c) correspond to Lemma 3.1. The swaps moving from 
G to its image G' is unique, therefore the probability of this transformation 
(the jumping probability from G to G' G) is: 

Prob(G G') := P(G'|G) = \ ■ (4.1) 

and 

Prob(G GO := P(G'|G) = \ • (4-2) 

(More precisely these are the probabilities that these vertex sets will be 
checked against making a swap.) The probability of transforming G to 
G' (or vice versa) is time-independent and symmetric. Therefore P is a 
symmetric matrix, where the entries in the main diagonal are non-zero, but 
(probably) different values. As we discussed it earlier (see Theorem 2.3), our 
Markov chain is irreducible (the state space is connected), and it is clearly 
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aperiodic, since it is lazy. Therefore, as it well known, our Markov process 
is reversible with the uniform distribution as the globally stable stationary 
distribution. 

Our main result is the following: 

Theorem 4.1. The Markov process defined above is rapidly mixing on each 
bipartite half-regular degree sequence with a forbidden star and a forbidden 
1 - factor. 

This result supersedes the previously mentioned three results ([9, 5, 13]) 
(except that this does not care of the actual mixing time, it just proves that 
it is polynomial). We want to add, however, that the friendly path method, 
described in [13], was not intended to handle half-regular bipartite degree 
sequences only but all bipartite degree sequences. Therefore we think that 
that method should not be completely neglected. 

There are several different methods to prove fast convergence of a Markov 

chain, here we apply a specialized version of Sinclair's seminal multicommod- 
ity flow method ([18]), the so called simplified Sinclair's method, developed 
in [13]^: 

Simplified Sinclair's method: We fix a half-regular star-l-factor restricted 
bipartite graphical sequence problem d^^. Our bipartite degree sequence is 
d = (a, b) where the vector a contains the degrees in class U while b 
contains the degrees in class W. (So all elements in a are the same, except 
maybe a(,s).) Therefore if X G G, then X is a simple bipartite graph 
([/, W; E{X)) and E{X) does not contain any element from F. The edge set 
E{G) corresponds to the possible swap operations. 

Sinclair's multicommodity flow method defines a bunch of paths (consec- 
utive sequences of swaps) for each realizations pair X and Y which transform 
realization X into Y. 

So consider two realizations X G G and F G G, and consider the sym- 
metric difference A = E{X)l\E{Y). In the bipartite graph 6 = ([/, VF; A) 
for each vertex v the number of adjacent X-edges (= E{X) \ EiY)) and 
the number of the adjacent F-edges are the same. Therefore, due to Euler 
classical reasoning, it can be decomposed into alternating circuits. 

The simplified Sinclair's multicommodity path method consists of two 
phases: In Phase 1 we decompose the symmetric difference A into alternat- 
ing circuits on all possible ways. In each cases we get an ordered sequence 
^1, W^2, ■ ■ ■ , Wk of circuits. (Usually there are a huge number of different 
decompositions.) Each circuit is endorsed with a fixed cyclic order. 

In Phase 2 each circuit Wi from the (ordered) decomposition derives 
one unique alternating cycles decomposition: Wi = C|, C2, . . . , C|,.. This 



^More precisely we need to slightly generalize it. It will be discussed in the Appendix. 
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decomposition is fully determined by the circuit and its well defined edge 
order. (Both construction algorithms are fully described in Section 5 of the 
paper [13], we do not discuss them here.) 

The ordered circuit decomposition together with the ordered cycle de- 
compositions of all circuits altogether provide a well defined ordered cycle 
decomposition Ci, . . . of A. 

This ordered cycle decomposition determines i—1 realizations Hi, . . . Hi-\ 
with the following property: if we use the notations Hq = X and = Y 
then for each j = 0, . . . ,f — 1 we have E{Hj)AE{Hj+i) = Cj+i. (It is im- 
portant to recognize that till this point we did not process even one swap 
operation! We just identified £ — 1 realizations which will be along our 
canonical path.) 

We will define a unique canonical path from X to Y determined by 
this circuit decomposition which uses these realizations Hj as milestones 
along the path. The canonical path will be X = Gq, . . . ,Gi, . . . , Cm = Y 
where each Gj can be derived from Gj-i with one valid swap operation, 
where we must have the following property: there are some increasing indices 
< ni < 722 < • • • < such that we have G„. = Hi. This, together the 
definitions of Hi means that 



The canonical path we are looking for has two important further properties: 
for each i < £ the constructed path Hi = Gq,Gi, . . . , G'^, = -ffj+i between 

Grii and G„-^-^ must satisfies that 

(G) m! <c-\Ci\ for a suitable constant c; 



(Q) for each j there is Kj G V{G) such that f Mx + My - Mq' , Mk^ ) < 



where the notations Mq stands for the usual bipartite adjacency matrix 

of G (this will be defined in details at the beginning of the next section), 
and d stands for the Hamming distance of two matrices of the same form, 
finally is a small constant. 

The current value of the auxiliary matrix Mx + My — Mqi together with 
the symmetric difference A, furthermore a small (polynomial) size parameter 
set, finally the vertices in G on which the canonical path under investiga- 
tion goes through uniquely determine the vertices X, Y and the path itself. 
Therefore it can be used to control certain features of the canonical path 
system. If the overall number of these auxiliary matrices are small (their 
number is smaller than a small polynomial of n multiplied with the number 
of possible realizations - as it is ensured by (^i)), then - as it was proved in 
[13] - our Markov chain is rapidly mixing. 




^2, 
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So in Phase 2 we have to build up our swap sequence between Hi and 
-fTj+i for all values i taking care for conditions (©) and (il). This will happen 
in the next Section. 

5. The construction of swap sequences between consecutive "mile- 
stones" 

Now we are going to implement our plan described above. At first we intro- 
duce some shorthand. Instead of Hi and -ff^+i we will use the names G and 
G'. These two graphs have almost the same edge set. More precisely 

{E{G) \ id n E{X))) U id n E{Y)) = E{G') 
{E{G') \ {Ci n E{Y))) U {Ci n E{X)) = E{G). 

Of course E{G)AE{G') = Ci also holds. We refer for the elements of Ci n 
E{X) as X-edges while the others are the F-edges. We denote the cycle 
itself as C, it has 2£ edges and its vertices are ui,wi,U2,W2,---, ug,wg. Since 
C has at least four vertices, therefore we may assume that ui ^ s (so ui is 
not the center of the forbidden star). Finally w.l.o.g. we may assume that 
the chord uiwi is an Y-edge (and, of course, w^ui is an X-edge). 

We are going to construct one by one the realizations G'j. We build 
our canonical path from C toward C and at any particular point the last 
constructed realization is denoted by Z. (At the beginning of the process we 
have Z = G.) We are looking for the next realization, denoted by Z'. 

Before we continue the discussion of the canonical path system, we have 
to introduce our control mechanism, mentioned in condition (ft). This aux- 
iliary structure originally was introduced by Kannan, Tetali and Vempala 
in [9]: 

For any particular realization G from G the matrix Mq denotes the 
adjacency matrix of the bipartite realization G where the columns and rows 
are indexed by the vertices of U and W resp. (Therefore the column sums 
are the same in each realization, except perhaps at column s.) Our indexing 
method is a bit unusual: the columns are numbered from left to right while 
the rows are numbered from bottom to the top. (Like in the Cartesian 
coordinate system.) This matrix is not necessarily symmetric, and elements 
Mi^i can be different from 0. 

For example if we consider the submatrix in Mq spanned hy ui, . . . ,U£ 
and wi, . . . ,wg then we have Mcii, i) = for i = while Mdi, i — 

1) = 1 (for i = 2, ...,£) and Mg(1,^) = 1. (So the first value gives the 
column, the second one gives the row.) The non-chords between vertices in 
the same vertex class are not considered at all, while non-chords which are 
forbidden are denoted by ^. As it is clear from the previous sentence, we 
will identify each chord or non-chord with the corresponding position in the 
matrix. 
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Our auxiliary structure is the matrix 

M{X + Y - Z) = Mx+My- Mz. 

By definition, each entry of a bipartite adjacency matrix is or 1 (or ^). 
Therefore only —1, 0, 1, 2 can be the "meaningful" entries of M. An entry is 
— 1 if the edge is missing from both X and Y but it exists in Z. This is 2 if 
the edge is missing from Z but exists in both X and Y. It is 1 if the edge 
exists in all three graphs {X, Y, Z) or it is there only in one of X and Y but 
not in Z. Finally it is if the edge is missing from all three graphs, or the 
edge exists in exactly one of X and Y and in Z. (Therefore if an edge exists 
in exactly one of X and Y then the corresponding chord in M is always 
or 1.) One more important, but easy fact is the following: 

Observation 5.1. The row and column sums of AI{X + Y — Z) are the 

same as row and column sums in Mx (or My or Mz)- O 

Next we will determine the swap sequence between G and G' through an it- 
erative algorithm. At the first iteration we check, step by step, the positions 
(ui, {ui.,w^), . . . , {ui,Wi) and take the smallest j for which {ui,Wi) is 
an actual edge in G. Since {ui,W(:) is an edge, therefore such i always exists. 
So we may face to the following configuration: 




edge chord, non-edge non-chord •'•.••■■unknown 

Figure 1; Sweeping a cycle 

We will call this {ui,Wi) chord as start-chord of our current sub-process 
and uiwi is the end-chord. We will sweep the alternating chords along 
the cycle from the start-edge WiUi (non-edge), UiWi-i (an edge) toward the 
end-edge wiui (non-edge) - switching their status in twos and fours. We 
check positions uiWi-i.,uiWi-2 (all are non-edges) and choose the first chord 
among them, we will call it the current-chord. (Since ui ^ s therefore we 
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never have to check more than two edges to find the first chord, and we need 
only one times to check two, since there is at most one non-chord adjacent 

to Ui.) 

Case 1: As we just explained the typical situation is that the current- 
chord is the "next" one, so when we start this is typically uiWi^i. As- 
sume that this is a chord. Then we can proceed with the swap operation 
Wi-iUi,WiUi uiWi-i,UiWi. We just produced the first "new" realization 
in our sequence, this is G'l. For the next swap operation this will be our new 
current realization. This operation will be called a single-step. 

In a realization Z we will call a chord bad, if its current status (being 
edge or non-edge) is different from its status in G (or, what is the same, 
in G', since they differ only on the chords along the cycle C). After the 
previous swap, we have two bad chords in G'^, namely uiWi-i and WiUi. 

Consider now the auxiliary matrix M{X + Y — Z) (here Z = G']). As 
we saw earlier, for each position outside the chords in C the status of that 
particular position in Z is the same as in X or F or in both. Accordingly, 
the corresponding matrix value is or 1. We call a position bad in M if this 
value is —1 or 2. (A bad position in M always corresponds to a bad chord.) 
Since in Case 1 we switch the start-chord into non-edge, it may become 2 
in M. (In case if in both X and Y it is an edge. Otherwise it is or 1, so in 
that case it is not a bad position.) The current-chord turned into an edge. 
If it is non-edge in both X and Y then the value becomes —1, otherwise it 
does not become a bad position. After this single-step, we have at most two 
bad positions in the matrix, at most one position with 2-value and at most 
one with —1- value. 

Case 2: If the position "below" the start-chord is a non-chord, then we 
cannot produce the previous swap. Then, however, the non-cdgc uiWi-2 is 
the current-chord. For sake of simplicity wc assume that i — 2 = 2 so we are 
in Figure 1. Consider now the alternating Cq cycle: ui,W2-,u^-,w^,U4^,Wi. It 
has altogether three vertex pairs which may be chords. We know already 
that u\W2, is a non-chord. If none of the three is chord, then this is an 
J^-compatiblc circular Cg-swap - and accordingly to the definitions we can 
swap it in one step. Again, we found the valid swap W2U^,w^Ui,WiUi =^ 
uiW2,U3Ws,U4U!4. After that we again have 2 bad chords, namely uiW2 and 
W4U1, and together we have at most two bad positions in the new M{X + 
Y — Z) with at most one 2-value and at most one — 1-value. 

Finally if one position, say ^21*4, is a chord then we can process this Cq 
with two swap operations. If this chord is, say, an actual edge, then we swap 
W2U4,,W4Ui uiW2,U4W4. After this we can take care for the W2,U3,W3,U4 
cycle. Along this sequence we never create more, than 3 bad chords: the first 
swap makes chords W2U4,W4Ui and U1W2 bad ones, and the second "cures" 
W2U4 but does not touch ui'W2 and W4U1. So along this swap sequence we 
have 3 bad chords, at the end we have only 2. On the other hand, if the chord 
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W2U4, is not an edge, then we can swap W2U^,W3U4 U3W3,U4W2, creating 
one bad edge, then taking care the four cycle ui,W2,U4,W4 we "cure" 'W2U4 
but we switch U1W2 and W4U1 into bad chords. We finished our double-step 
along the cycle. 

In a double-step we make at most three bad chords. When the first swap 
uses three chords along the cycle then we may have at most one bad chord 
(with M- value or —1) and then the next swap switches back the chord into 
its original status, and makes two new bad chords (with at most one 2-value 
and one — 1-value). When the first swap uses only one chord from the cycle, 
then it makes three bad chords (changing two chords into non-edge and one 
into edge), therefore it may make at most two 2- values and one —1-value. 
After the second swap there will be only two bad chords, with at most one 
2-value, and at most one —1-value. 

When only the third position corresponds to a chord in our Cq then after 
the first swap we may have two — 1-values and one 2-value. However, again 
after the next swap we will have at most one of both types. 

Remark 5.2. When two realizations are one swap apart (so they are ad- 
jacent in G) then we say that their auxiliary matrices are at swap-distance 
one. Since one swap changes four positions of the matrix, therefore the 
Hamming distance of these matrices is 4. 

Finishing our single- or double-step the previous current-chord becomes the 
new start-chord and we look for the new current-chord. Then we repeat our 
procedure. There is only one important point to be mentioned: along the 
step, the start-chord switches back into its original status, so it will not be 
a bad chord anymore. So even if we face a double-step the number of bad 
chords never will be bigger than three (together with the chord WiUi which 
is still in the wrong status, so it is bad), and we have always at most two 
2- values and at most one —1-value in M{X + Y — Z). 

When our current-chord becomes to W1U2 then the last step will switch 
back the last start-chord into its correct status, and the last current-chord 
cannot be in bad status. So, when we finish our sweep from uiWi to wiui 
at the end we will have only one bad chord (with a possible 2-value in M). 
This concludes the first iteration of our algorithm. 

For the next iteration we seeks a new start-chord between WiUi and W£Ui 
and chord WiUi becomes the new end-chord. We will repeat our sweeping 
process for this setup, and we will repeat it as long as all chords will be 
processed, so we fond the entire realization sequence from G to G' . If in the 
first sweep we had a double-step, then it will never occur later, so altogether 
with the bad (new) end-chord we never have more than three bad chords, 
with at most two 2- values and at most one —1-value. 

However, if the double-step occurs sometimes later, for example in the 
second sweep, then we face to the following situation: if we perform a circular 
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Cg-swap, then there cannot be any problem. So we may assume that there 
is a chord in our Cg, suitable for a swap. If this chord is a non-edge, then the 
swap around it produces one bad chord, and at most one bad position in M. 
The only remaining case when that chord is an edge. After the first swap 
there will be four bad chords, and there may be at most three 2-valucs and 
at most one —1 value. However after the next swap (finishing the double 
step) we annihilate one of the 2- values, and after that swap there are at most 
two 2- values and at most — 1-value along the entire swap sequence. When 
we finish our second sweep, then chord WiUi will be switched back into its 
original status, it will not be bad anymore. 

We apply iteratively the same algorithm, and after at most £ sweep 
sequence, we will process the entire cycle C. This finishes the construction 
of the required swap sequence (and the required realization sequence). □ 

Meanwhile we also proved the following important observation: 

Lemma 5.3. Along our procedure each occurring auxiliary matrix M(X + 
Y — Z) is at most swap- distance one from a matrix with at most three bad 
positions: with at most two 2-values and with at most one —1-value in the 
same column, which does not coincide with the center of the forbidden star. 

6. The analysis of the swap sequences between "milestones" 

What remains is to show that the defined swap sequences between Hi and 
Hi^i satisfy the properties (O) and (ft) of the simplified Sinclair's method. 
The first one is easier to see, since we can process a cycle of length 2£ in 
£ — 1 swaps. Therefore the derived constant c in (G) is actually 1. 

We introduce the switch operation on 0/1 matrices with forbidden po- 
sitions: we fix the four corners of a submatrix (none of them is forbidden), 
and we add 1 to two corners in a diagonal, and add —1 to the corners on 
the other diagonal. This operation clearly does not change the column and 
row sums of the matrix. For example if we consider the matrix Mq of a 
realization of our d"^ and make a valid swap operation, than it looks like as 
a switch in this matrix. The next statement is trivial but very useful: 

Lemma 6.1. If two matrices have switch-distance 1, then their Hamming 
distance is 4. Consequently if the switch- distance is c then the Hamming 
distance is bounded by 4c. 

We will prove now that property (fi) holds for our auxiliary matrices: 

Theorem 6.2. For any realizations X and Y furthermore for any realiza- 
tion Z on a swap sequence from X to Y there exists a realization K such 
that 
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Due to Lemmas 5.3 and 6.1 it is enough to show that: 

Lemma 6.3. Any matrix M{X + Y — Z) with constant column sums (this 
not necessarily holds for the center of the forbidden star) and at most three 

bad positions (where there are at most two 2-values and at most one — 1- 
value) can he transformed into a valid adjacency matrix with at most 
three switch operations. 

Proof. Consider now a certain M which is not necessarily a vahd adjacency 
matrix of a reahzation. We wih show pictures about the submatrix in this 
matrix which describes the current alternating cycle C. We choose a subma- 
trix, where the center s of the forbidden star is in the first column. (We 
choose this submatrix as an illustration tool, but we still consider the en- 
tire matrix to work with.) We know that this matrix contains at most two 
2-values and at most one 1-value. All these positions are adjacent to the 
center u\ of our sweeping sequence (see Figure 1), so they are in the same 
column. 

For simplicity from now on we will denote the center of the sweep as well 
the column u. The forbidden positions are denoted with Any column 
(except column 1) may contain at most one of them, and any row may 
contain at most two of them. Finally in our pictures the character o stands 
for a character which we are not interested in. That is, it can be or 1 or 

We will distinguish cases, depending on the occurring of values 2 and 

-1. 

Case 1. Column u has one bad position, which can be —1 or 2, or it has 
two 2-values. Consider at first the subcase when M[«i(;] = —1. By definition 
that means that chord uw is an edge in Z but non-edge in both X and Y. 
So vertex w W has at least one adjacent edge, therefore the row-sum in 
its row is at least 1. Therefore there are at least two positions in row w 
with entries 1. They are in column ui and U2. At least one of them, say ui, 
differs from s. Since the column sums are constant, therefore there exists at 
least two rows wi such that M[utt;i] = 1 while M[uiii;i] = or However, 
there can be at most one forbidden position in ui, so in at least one of the 
rows, the entry is 0. Using these positions for the corresponding switch it 
eliminates the bad position without creating a new one. (See Figure 2.) 
Before we continue, we prove an important observation: 
Observation If w belongs to the alternating cycle C and M['ut«] = 2 then 
row w contains at least two 0-values. 

Indeed, there are a forbidden chords in row w. Since w is in an alternating 
cycle, therefore d{w) < |?7| — a — 1. Therefore the sum of row w in M{X + 
Y — Z) < \U\ — a — l. But it contains a 2 and it does not contain -1 therefore 
there are at least two O's in it. □ 
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Figure 2: Case 1 "J = forbidden o = 0/l/>J< 

When the single bad vahic in M is 2 then, due to our previous Obser- 
vation, in its row there are two O's. And with them one can repeat the 
reasoning which we used about the unique —1- value. 

Finally, when there are two 2- values which raises a very similar situation. 
Here we can do the same procedure independently on both rows. In this 
case, however, we need two switch operations. 

Case 2. Here we assume that there is one 2- value and one — 1-value in 
column u. For example M[uw\\ = 2 and M[uw2\ = —1. Again, in row W2 
there are at least two 1-values. 

Case 2a Assume at first that we have ui £ U s.t. M[uiW2] = 1 and 
M[uwi] ^b. Then the corresponding switch will produce M\u,\W\\ = 1/2 
while the three positions are or 1. (See Figure 3.) If now M\u,\W\\ = 2 
then we are back to Case 1, and one more switch eliminates the last bad 
position as well. So we needed at most two switches. 
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Figure 3: Case 2a 5< = forbidden 



:0/l/^< 



Case 2b It can happen, that there are only two 1-values in row W2 and both 
are facing with forbidden positions in row wi. Then at least one in row W2 
faces a chord in row W2- (See Figure 4) The appropriate switch kills 2 bad 
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Figure 4: Case 2b ^ = forbidden o = 0/1/^ 



chords and can make at most one —1- value. We are ready or we are back 
to Case 1. 

Case 3. Finally suppose that there are three bad positions, two 2-values at 
positions uwi and u'W2 and one — 1-value at position uws. Now both rows 
wi and W2 contains at least two O's. If any of them faces a 1 in row W3 then 
an appropriate switch annihilates one 2 and one —1 and does not create new 
bad position. We are back to Case 1. Altogether we need two switches. 

If this is not the case then we consider the following: assume that 
MfuifWi] = 0. Since the column sums are the same, and we assumed that 
M[uiW3] = therefore there exists a row W4 s.t. M[uiW4] = 1 while 
M['ui(;4] = 0. Then we can switch off this 2-value without making a new 
bad position. After that we are back to Case 2. Altogether this requires at 
most three switches. We finished the proof of Lemma 6.3. 

If this is not the case then we consider the following: assume that 
M[uiwi] = 0. The column sums are the same, and we assumed that M[uiW3] 
= or ^. Therefore the difference between column sums in u and ui is 1 due 
to rows wi and W3, and the difference increase at least 1 for row W2, where 
against a 2-value in column u there is either 1 or in column ui. Therefore 
there exists at least two further rows, where there is a 1 in column ui against 
a or ^ in column u. Since column u can contain at most one one of the 
rows must contain a 0. Let it be denoted by W4. Hence MluiWi] = 1 while 
M[uw4] = 0. Then we can switch off this 2-value without making a new bad 
position. After that we are back to Case 2. Altogether this requires at most 
three switches. We finished the proof of Lemma 6.3. □ 

In turn this proves Theorem 6.2, so our Markov chain is rapidly mixing 
as Theorem 4.1 stated. 
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7. Self-reduced counting problem 

A decision problem is in NP if a non-deterministic Turing Machine can 
solve it in polynomial time. An equivalent definition is that there exists a 
witness proving the yes answer to the question which witness can be verified 
in polynomial time. A counting problem is in #P if it asks for the number 
of those witnesses of a problem from NP that can be verified in polynomial 
time (it might happen that not all witnesses are verifiable in polynomial 
time) . 

Two complexity classes, FPRAS and FPAUS, concern the approximabil- 
ity of counting problems. Here we give only narrative descriptions of these 
complexity classes, the detailed definitions can be found, for example, in [8]. 

A counting problem from #P is in FPRAS (FtiIIv Polynomial Random- 
ized Approximation Scheme) if the number of solutions can be quickly es- 
timated with a randomized algorithm such that the estimation has a small 
relative error with very high probability. 

A counting problem from ^^P is in FPAUS (Fully Polynomial Almost 
Uniform Sampler) if the solutions can be sampled quickly with a randomized 
algorithm that generates samples following a distribution being very close 
to the uniform one. 

It is easy to see that a counting problem is in FPAUS if there is a rapidly 
mixing Markov chain for which 

• a starting state can be generated in polynomial running time; 

• one step in the Markov chain can be conducted in polynomial running 
time; and 

• the relaxation time of the Markov chain grows only polynomially with 
the size of the problem. 

The Markov chain we gave the star-l-factor problem satisfies all these re- 
quirements. 

Jerrum, Valiant and Vazirani proved that any self-reducible counting 
problem is in FPRAS iff it is in FPAUS [8]. A counting problem is self- 
reducible if the solutions for any problem instance can be generated re- 
cursively such that after each step in the recursion, the remaining task is 
another problem instance from the same problem, and the number of possi- 
ble branches at each recursion step is polynomially bounded by the size of 
the problem instance. 

Clearly, a graph with prescribed degree sequence can be built recursively 
by telling the neighbors of a node at each step, then removing the node 
in question and reducing the degrees of the selected neighbors. However, 
this type of recursion does not satisfy all the requirement for being self- 
redTiciblc since there might be exponentially many possibilities how to select 
the neighbors of a given vertex. 
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On the other hand, the degree sequence problem with a forbidden one 
factor and star tree is a self-reducible counting problem. Indeed, consider 
the center of the (possibly empty) star, s € U, and the vertex v & V with 
the smallest index for which {s,v) is a chord. Any solution for the current 
problem instance belongs to one of the following two cases: 

• The chord {s,v) is not present in the solution. In that case, extend 
the size of the star by adding chord {s,v) to the forbidden set, and 
do not change the degrees. This is another problem instance from 
the star+factor problem, whose solutions are the continuations of the 
original problem belonging to this case. 

• The chord (s, v) is present in the solution. In that case, extend the size 
of the star by adding chord (s, v) to the forbidden set, and decrease 
both ds and by one. The new degree sequence is still a half-regular, 

bipartite star+factor restricted degree sequence, and the solutions of 
this new problem extended with the previously decided step provide 
solutions of the original problem. 

Since the star+factor counting problem is a self reducible counting problem, 
it is in FPRAS as it is in FPAUS. 
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Appendix: the simplified Sinclair's method 

In this paper the fast mixing nature of our MCMC method was proved 
through the application of the simplified Sinclair's method, developed in 
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[13]. To do so properly it requires a slight generalization of the original 
method. 

The method takes two realizations X and Y of the same degree sequence. 
It considers all possible ordered circuit decompositions of the symmetric dif- 
ference of the edge sets, then it uniquely decomposes each such decomposi- 
tion into an ordered sequence C = Ci, . . . , Cm of oriented cycles. Based on 
this latter decomposition the method determines a well defined unique path 
between X and Y in the Markov chain G. 

For that end the method defines first a sequence of "milestones". These 
are different realizations X = Hq, Hi, . . . , Hm-i, Hm = Y of the degree 
sequence where the edge set of any two consecutive realizations Hi-i,Hi 
differ exactly in the edges along the cycle Q. (Until this point no swap 
operation happened.) 

In the next phase for any particular i = 0, . . . ,m — I the method de- 
termines a sequence of valid swap operations transforming Hi-i into - 
describing a unique path Zq, Zi, . . . between -ffj-i and in the Markov 
chain G. This sequence of course heavily depends on the available swap oper- 
ations. In paper [13] these are the usual (bipartite) swap operations. In the 
current paper these are the restricted swap operations. These operations, 
while exchanging chords in the realizations along the alternating cycle Cj, 
also use some further chords. Therefore the edge set of any Zi is not com- 
pletely contained by E{X) U E{Y), there exist a small number of edges in 
Zi which are non-edges in X and in y, or non-edges in Zi but edges in X 
and Y. If Zi is between the milestones H^n and Hm+i, then Cj for j ^ m 
alternates in Zi, and Q alternates with a "small error": there is a very small 
number of vertices where the alternation does not hold. 

Along the process the simplified Sinclair's method requires (see the paper 
[13], Section 5, (F)(c) ) that this number must be small. In the original 
application this number is actually one. Here, as we saw in Section 5, this 
number is three: that many bad chords may occur after any particular RSO. 
As we saw all these chords are adjacent to the same vertex ui. 

These numbers are used by the method to determine the size of a pa- 
rameter set B. This parameters set must have a polynomial size. When we 
have one bad chord, then it is determined by its end points - there are at 
most possibilities for them. This provides an multiplicative factor to 
the size of B. When we have at most three bad chords, then they can be 
chosen at most ways: point ui is fixed {n different choices), while the 
other three end points can be chosen at most independent ways. Alto- 
gether it provides an at most multiplicative factor to the size of B. This 
remark finishes the proof of the simplified Sinclair's method for the case of 
these restricted swap operations. 
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